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ABSTRACT 

Gamma-ray emission from the Crab Nebula has been recently shown to be unsteady. In this paper, we study the 
flux and spectral variability of the Crab above 100 MeV on different timescales ranging from days to weeks. In 
addition to the four main intense and day-long flares detected by AGILE and Fermi-LAT between Sept. 2007 and 
Sept. 2012, we find evidence for week- long and less intense episodes of enhanced gamma-ray emission that we call 
"waves" . Statistically significant "waves" show timescales of 1-2 weeks, and can occur by themselves or in association 
with shorter flares. We present a refined flux and spectral analysis of the Sept. - Oct. 2007 gamma-ray enhancement 
episode detected by AGILE that shows both "wave" and flaring behavior. We extend our analysis to the publicly 
available Fermi-LAT dataset and show that several additional "wave" episodes can be identifled. We discuss the 
spectral properties of the September 2007 "wave" /flare event and show that the physical properties of the "waves" 
are intermediate between steady and flaring states. Plasma instabilities inducing "waves" appear to involve spatial 
distances I ~ 10^^ cm and enhanced magnetic fields B ~ (0.5 — 1) mC. Day-long flares are characterized by smaller 
distances and larger local magnetic fields. Typically, the deduced total energy associated with the "wave" phenomenon 
{Eyj ~ 10''^ erg, where is the kinetic energy of the emitting particles) is comparable with that associated to the 
flares, and can reach a few percent of the total available pulsar spindown energy. Most likely, flares and waves are the 
product of the same class of plasma instabilities that we show acting on different timescales and radiation intensities. 
Subject headings: 



1. INTRODUCTION 

The Crab Nebula (the remnant of a Supernova explo- 
sion witnessed by Chinese astronomers in 1054) is pow- 
ered by a very powerful pulsar (of period P = 0.33 ms, 
and spindown luminosity Lsd ~ 5x 10'^* ergs^^) (see e.g., 
Hester 2008). The pulsar is energizing the whole system 
through the interaction of the particle and wave output 
within the surrounding Nebula (of average magnetic flcld 

200/zG). The resulting unpulsed emission from radio 
to gamma rays up to 100 MeV is interpreted as syn- 
chrotron radiation from at least two populations of elec- 
trons / positrons energized by the pulsar wind and by sur- 
rounding shocks or plasma instabilities (e.g., Atoyan & 
Aharonian 1996, Meyer et al. 2010). The optical and X- 
ray brightness enhancements observed in the inner Neb- 
ula, known as "wisps" , "knots" , and the "anvil" aligned 
with the pulsar "jet" (Scargle 1969; Hester 1995, 2008; 
Weisskopf 2000), show flux variations on timescales of 
weeks or months. On the other hand, the average un- 
pulsed emission from the Crab Nebula was always con- 
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sidered essentially stable. The surprising discovery 
by the AGILE satellite of variable gamma-ray emission 
from the Crab Nebu la in Sept. 2010 (Tftyani et al. 2010l 
Tavani et al. 2011a), and the Fermi-LAT confirmation 



(I'Buehler et al. 20101 lAbdo et al.^OTTI) started a new 
era of investigation of the Crab system. As of Sept. 2012 
we know of four major gamma-ray flares from the Crab 
Nebula detected by the AGILE Gamma- Ray Imaging De- 
tector (GRID) and Fermi-LAT: (1) the Sept-Oct. 2007 
event, (2) the Feb. 2009 event, (3) the Sept. 2010, and 
(4) the "super-flare" event of Apr. 2011 (Buehler et al. 
2011; Tavani et al. 2011b; Hays et al. 2011; Striani et 
al. 2011a; Striani et al. 2011b; Buehler et al. 2012). 

In this paper we address the issue of the gamma-ray 
variability of the Crab Nebula on different timescales 
ranging from days to weeks. We then enlarge the pa- 
rameter space sampled by previous investigations espe- 
cially for the search of statistically signiflcant enhanced 
emission on timescales of 1-2 weeks. Sect. 2 presents a 
brief overview of the current knowledge on Crab's main 
gamma-ray flares. Sect. 3 presents the results of a search 
of 7-ray enhanced emission on timescales of weeks in the 
AGILE database. We also discuss in detail the Sept. -Oct. 
2007 event detected by AGILE which shows a strong ev- 
idence of short timescale flaring as well as substantial 
emission on longer timescales of order of 1 week. Sect. 
4 presents the results of our search for long and short 
timescale enhanced emission in the available Fermi-LAT 
data. In both the 7-ray telescopes data we find strong 
evidence of week-long 7-ray enhanced emission episodes 
at intermediate peak intensities that we call "waves". 
Sect. 5 presents the physical implications of our findings 
in terms of a synchrotron emission model. We discuss in 
Sect. 6 the main implications of our work. 
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TABLE 1 

Table of the flares (F > 700 x 10~* phcm~^ s"-*^) of the Crab Nebula found in the AGILE and Fermi data from Sept. 2007. 





Name 


MJD 


n (hr) 


T2 (hr) 


Peak Flux 


B{mG) 


7* (10^) 


I (10^^ cm) 


9007 


Pi 
^1 


O^OOi .0 


99 -1-11 


1 n -t- ^ 


1 nnn -i- i f^n 
iuuu m ±ou 


in 9 


Z.D — ^.o 


19 


(AGILE) 


F2 


54382.5 


14 ±7 


6±3 


1400 ± 200 


1.1 - 2.1 


2.3 - 4.3 


0.8 - 2.2 




F3 


54383.7 


11 ± 5 


14 ± 7 


900 ± 150 


1.0 - 2.0 


2.6-4.8 


0.8 - 1.7 


2009 


F4 


54865.8 


10 ±5 


20± 10 


700 ± 140 


0.7- 1.3 


2.6-4.8 


0.6 - 1.6 


(FERMI) 


Fs 


54869.2 


10 ±5 


22 ± 11 


830 ± 90 


0.8 - 1.4 


2.6-4.8 


0.6 - 1.6 


2010 


Fa 


55457.8 


8±4 


22 ± 11 


850 ± 130 


0.7- 1.3 


2.5-4.7 


0.5 - 1.3 


(AGILE & 


Fv 


55459.8 


6±3 


6±3 


1000 ± 100 


1.4 - 2.6 


2.6 - 4.8 


0.3 - 0.9 


FERMI) 


Fs 


55461.9 


19± 10 


8±4 


750 ± 110 


0.8 - 1.4 


2.5 - 4.8 


0.9 - 3.1 


2011 


F9 


55665.0 


9±5 


9±5 


1480 ± 80 


1.2 - 2.2 


2.8 - 5.0 


0.5 - 1.5 


(FERMI & 


Fio 


55667.3 


10 ±5 


24± 12 


2200 ± 85 


1.3 - 2.3 


2.7 - 4.9 


0.6 - 1.6 


AGILE) 



















The timescales ri and T2 are the rise and decay timescales of the flares modelled with an exponential fit respectively. The characteristic 
length of the emitting region is deduced from the relation I = cSti. The Lorentz factor 7* characterizes the adopted model of the 
accelerated particle distribution function dn/d'y = K/a ■ 6(7 — 7*), where K is defined in the spherical approximation, a = 1 in the 
spherical case, and a < 1 for cylindrical or pancake-like volumes reproducing the current sheet geometry. The peak photon flux above 
100 MeV is measured in units of 10~* phcm~^ s~^. These parameters are obtained, by means of a multi-parameter flt, from the following 
quantities (in the observer frame): the position of the peak photon energy, Ep oc S^*^B, the peak emitted power vF oc &'^K / al^ B^'y*^ , 
the rise time ri = l/{cS), and the cooling time T2 = 8.9 X 10^/[(iJ/Gauss)^ 7*^5]i assuming (5 = 1 (see text). 




Fig. 1. — Gamma-ray lightcurves above 100 MeV (12-hr time bins) from the Crab (pulsar plus Nebula) detected by AGILE and Fermi- 
LAT. From top to bottom, the Sept. - Oct. 2007 event (AGILE data), the Feb. 2009 event (Fermi-LAT data), the Sept. 2010 event 
(Fermi-LAT data), and the Apr. 2011 event (Fermi-LAT data). 
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Fig. 2. — Gamma-ray 1-day binned lightcurve above 100 MeV of the five AGILE observations of the Crab Nebula in pointing mode, 
from Sept. 2007 to Oct. 2009. The gamma-ray fiux enhancement during the September-October 2007 pointing (MJD = 54366 - 54386) 
is evident. 
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Fig. 3. — Plot of the values (each calculated for 7-day time intervals based on 1-day binned data) as a function of time for the AGILE 
data on the Crab covering the period Sept. 2007/Oct. 2009 in pointing mode. The red dashed lines indicate the corresponding to the 
3(T, 4a and 5cr confidence level. 
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Fig. 4. — Lightcurve (1-day bin) of the Sept.- Oct. 2007 Crab Nebula flare detected by AGILE. In the inset the 12-hr bin hghtcurve 
around the flare. This episode is characterized by a very strong variability, with waves (black line, marked with a W) and flares (red line, 
marked with an F.) 



2. OVERVIEW OF THE MAIN GAMMA-RAY FLARES 

Four major episodes of intense gamma-ray flaring from 
the Crab Ne bula fiave been dete cte d by AGILE and 
Fermi-LAT (ITavani et al. 201 lal lAbdo et al. 20111 
[Striani et al. 2011b; IVittorini et al. 20111 

Euehlcr ct al. 2012). The definition of a "flare" 
adopted in this paper is that of a single gamma-ray 
enhancement event with a risetime ti < 1 day and 
flux F > 700 X lO^^phcm-^s-i above 100 MeV. 
Table 1 summarizes the flaring events that we find 
by considering the AGILE and the publicly available 
Fermi-LAT database. These events show a complex 
time structure, being composed of several sub-flares 
that we classify as Fn', with n' a progressive number. 
Fig. 1 summarizes the four major flaring episodes with 
the same temporal and flux scales. The colored curves 
are indicative of the flaring behavior that in most cases 
can be represented by an exponential fit^^(e.g., Norris 
et al. 1996), characterized by rising (n) and decay {T2) 
timescales as given in Table 1. We also report in Table 
1 the data fitting physical parameters (average local 
magnetic field B, typical particle Lorentz factor 7*, and 
characteristic size / of the emitting region of the flaring 
episodes). We determine these parameters and their 
uncertainties from the time constants ti and T2 (see 
discussion below). Whenever applicable (2010 and 2011 
events), the AGILE and Fermi-LAT data are consistent 
both in flux and spectral properties. 

Motivated by the claim of fast variability in the 
September 2010 event (Balbo et al. 2011) and the dra- 
matic detection of the short timescale variability in the 
Crab gamma-ray flare of April 2011 (Buehler et al. 2011), 



We use the fltting function 

r F6-hAexp(-| 
\F6-fAexp(-| 



tp\/rji) for t < tp 
■ *pl/o"2) for t > tp 



where Ff, is the flux baseline, A a flux parameter, tp is the peak 
time, (Ti and 172 the rise and decay time constants. The rise and 
decay times, half to maximum amplitude, are obtained as ti 2 = 
[ln(2)](7i,2. 



we revisited the analysis of the Crab flare detected by 
Fermi in Feb. 2009, and the AGILE and Fermi analy- 
sis of the September 2010 event (12-hr bin lightcurve for 
E > 100 MeV, first two panels of Fig. fl|. This revised 
analysis shows for both the 2009 and 2010 events a se- 
quence of three flares, that could not be previously appre- 
ciated with a 2-day time bin. Regarding the September 
2010 event, both the AGILE data and the Fermi-LAT 
data agree quite well, and are consistent with the analy- 
sis presented in Balbo et al. 2011. 

3. THE SEPTEMBER-OCTOBER 2007 EVENT DETECTED 
BY AGILE 

The Crab pulsar plus Nebula is a primary source for 
gamma-ray calibration, and AGILE pointed at the source 
several times during the pointing mode phase from July 
2007 until Oct. 2009 for different geometries and off- 
axis angles. The ideal periods during the year for AG- 
ILE pointings of the Crab region are September-October 
and March-April as determined by the solar panel con- 
straints. Fig. [2] shows the overall gamma-ray lightcurve 
collecting all the available AGILE observations which in- 
cludes several pointings in 2007, 2008 and 2009. Typi- 
cally, pointing observations lasted ^10 days, except for 
the initial pointing in September-October 2007 that was 
substantially longer. 

During this first pointing, indeed, we detected an 
episode of enhanced gamma-ray emission. A prelim- 
inary lightcurve of this episode, performed with a 1- 
day bin, was presented in Tavani et al. 2011a. In 
that study we reported a quite long episode (about 2 
weeks) of enhanced emission with a 1-day peak flux of 
Fp = (890 ± 110) X 10~** phcm-2 g-i above 100 MeV. 

Fig. [2| shows the 1-day hghtcurve {E > lOOMeV) of 
the Crab (pulsar plus Nebula) extended to the five ob- 
servations of the Crab Nebula performed by AGILE dur- 
ing the pointing mode. In order to search for sta- 
tistically significant enhanced gamma-ray emission on 
timescales of weeks, we calculated from the data of Fig. [2] 
the 7-day curve^'* (Fig. [3| in the null hypothesis 

^'^ For the AGILE data, we started with the 1-day binned flux 
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that the Crab Nebula is constant at its average flux 
{Fs = 220 X lO^^phcm^^s^i in the Catalog of Pittori 
et al. 2009). 

We then calculated the probability of obtaining a given 
(considering the degrees of freedom) in the null hy- 
pothesis. If Xobs the observed value of x^, p(x^in) is 
the probability of obtaining a value of > Xobs, with 
n degrees of freedom. For 7 degrees of freedom, the 
horizontal dashed lines in Fig. [3] indicate the value of 
X^ (quantiles) corresponding to the 3ct, 4cr and Scr confi- 
dence level, respectively. We see that the post-October, 
2007 four observations of the Crab Nebula performed 
by AGILE in the pointing mode are compatible with the 
average emission within la. However, three episodes 
during the first observation (from MJD ~ 54360 to MJD 
~ 54395) are above 5a with respect to the Crab Nebula 
average emission. 

The lightcurve extending for 20 days including the first 
pointing, from Sept. 24, 2007 to Oct. 13, 2007, is shown 
in Figure|4] In the inset we show a zoom of the lightcurve 
focused on the short variability episodes. Three events 
stand out above 5 sigma (for a 7-day time bin) in Fig. 3. 
The most significant corresponds to a flaring sequence, 
and two other events can be attributed to a less intense 
but long timescale enhanced emission anticipating the 
flares, with an average flux of ~ 450 x 10^* phcm~^ 
and a rise and decay time of the order of several days. 

These slow components of enhanced 7-ray emission 
show features different than those of flares that typically 
have rise and decay times of the order of 12-24 hr, and 
peak fluxes ranging from Fp^^ ~ 800 x 10~* phcm~^ s~^ 
up to Fp^io ~ 2500 X 10~^ phcm~^ s^^ (as for the Crab 
super-fiare of Apr. 2011). 

These episodes of slow enhanced emission are signif- 
icant (see discussion below): we call them waves. We 
indicate with Wi the emission from MJD ~ 54367 to 
MJD ~ 54374. From MJD ~ 54376 to MJD ~ 54382 
we observe another event of enhanced emission, that we 
interpret as a second "wave" with a small flare superim- 
posed at MJD ~ 54381. We name this second region W2. 
The 12-hr lightcurve in the inset of Fig. |4] shows that the 
2007 peak intensity event, that in our previous analy- 
sis ([Tavani et al. 2011aP appeared unresolved, is actually 
composed of three different flaring components, that we 
indicate by Fi, F2, and F3. Peak fluxes, rise times and 
decay times (estimated with an exponential flt) for Fi, 
F2, and F3 and for Wi and W2 are presented in Tab. [T] 
and Tab. [3l In order to take into account the number of 
1-day maps (trials) carried out in our search for enhanced 
emission, we also calculated the post-trial significance for 
Wi and W2, where the post-trial probability is given by 
Ppost = 1 — (1 — p)^' with p the pre-trial probability 
obtained from the x^ test, and Nt the number of trials. 
AGILE in pointing mode observed the Crab for a total 
time of ^ 50 days; considering that the value is cal- 
culated over 7 days, the number of trials is Nt ~ 7. We 
find that the post-trial significance is Ppost > 5a for both 
Wi and W2. 

The spectral properties of the Wi and W2 waves and 

data, and then calculated the values summing over 7 days, = 
^l—i (ZoC?i^iiii_ ^ where Fo{i) is the i-th observed flux, is the 
steady state flux, and Ci is the i-th flux error. 
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Fig. 5. — AGILE-GRID gamma-ray pulsar-subtracted spectrum 
of the Crab Nebula Wi (integrating from MJD = 54368 to MJD 
= 54374, top panel), W2 (integrating from MJD = 54376 to MJD 
= 54381, central panel) and the major flare F2 (12-hr integration, 
from MJD = 54382.5 to MJD = 54383.0, bottom panel). 



of the flares in September-October 2007 are interesting. 
Fig. [5] shows the spectral analysis of Wi (obtained inte- 
grating from MJD = 54368 to MJD = 54374), W2 (ob- 
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tained integrating from MJD = 54376 to MJD = 54381), 
and the spectral analysis of F2 obtained with a 12-hr in- 
tegration. (We notice that the spectrum of W2 might be 
contaminated by presence of the first flare, Fi). 

We find that the differential particle energy distribu- 
tion function (per unit volume) can be described by a 
monoenergetic function, dn/dj = K / a - — ^*), where 
7 is the particle Lorentz factor, and 7* is the monochro- 
matic value of the particle energy. The constant K is 
defined in the spherical approximation, with a = 1 in 
the spherical case, and a < 1 for cylindrical or pancake- 
like volumes. A power-law distribution and/or a rela- 
tivistic Maxwellian distribution were also shown to be 
consistent with the flaring data (Tavani et al., 2011a, 
Striani et al. 2011b). We adopt here a monochromatic 
distribution that deconvolved with the synchrotron emis- 
sivity leads to an emitted spectrum practically indistin- 
guishable from the relativistic Maxwellian shape (e.g., 
Striani et al. 2011b, Buehler et al. 2012). This dis- 
tribution is in agreement with all available gamma-ray 
data (for both flaring and "wave" behavior, see below), 
and reflects an important property of the flaring Crab 
acceleration process (Tavani 2013). In our model we 
have five free physical parameters: the Lorentz factor 
7*, the local magnetic field S, the electron density con- 
stant the dimension of the emitting region I, and 
the Doppler factor 5. The values of these parameters are 
obtained, by means of a multi-parameter fit, from the fol- 
lowing quantities (in the observer frame): the position of 
the peak photon energy, Ej, cx S'y*'^B, the peak emitted 
power lyF cx the rise time ti = l/{cS), 

and the cooling time T2 = 8.9 X 10V[(5/Gauss)2 7*5]. 
We fix the Doppler factor at the value 5 = 1 (relax- 
ing this condition leads to slightly different constraints 
for B and 7* that can be easily calculated without al- 
tering the main conclusions of our paper). We deter- 
mine the characteristic timescales ti and T2 by a 3- 
parameter model (see note 1), and from those values we 
deduce the other physical quantities. We find for Wi 
a magnetic field B = (0.8 ± 0.2) mG, a Lorentz factor 
7* = (4 ± 1) X 10^, and a typical emitting length range 
/ = (0.5-1.5) X lO^^ cm. For F2, we find B = (1.5±0.5) 
mG, 7* = (3 ± 1) X 10^ and I = (1.5 ± 0.7) x lO^^ cm. 
In our model, the total number of accelerated particles 
producing the gamma-ray wave/flaring behavior is in the 
range iV - (1 - 3) • 1038(Ari/47r), with An the sohd an- 
gle of the 7-ray emission. It is interesting to note that 
besides the differing values of the magnetic fleld and par- 
ticle densities, the typical Lorentz factor and total num- 
ber of radiating particles are similar for the "wave" Wi 
and the flare F2. We find that this is a typical behavior 
of the transient gamma-ray emission that appears to be 
well represented by a monochromatic particle distribu- 
tion function with 7* ~ (3 — 5) x 10^. 

4. SEARCH FOR ENHANCED GAMMA-RAY EMISSION IN 
THE FERMI-LAT DATA 

Motivated by the waves found in the AGILE data, we 
searched for a similar type of enhanced gamma-ray emis- 
sion in the publicly available Fermi-LAT data. Fig. [6] 
shows the lightcurve of the Crab (pulsar -|- Nebula) dur- 
ing the period Sept. 2008 - May 2012 obtained by a 
standard unbinned likelihood analysis of the Fermi data 
(2-day bin). The three major flares from the Crab Neb- 



ula detected by Fermi-LAT, and the new low-intensity 
event recently announced in Jul. 2012 in ATel #4239 
(Ojha et al. 2012) are recognizable at MJD - 54869, 
MJD ~ 55459, MJD - 55667 and MJD - 56113. 

Starting with the 2-day bin lightcurve we calculated 
the distribution based on 8-day integrations, and 
tested the null hypothesis for a source with a constant 
average flux Fp = (296 ±2.5) x 10~^ phcm~^ s~^ above 
100 MeV for the pulsar plus Nebula signal^*^. We used 
the same procedure employed for the study of the AG- 
ILE data^^. In Fig. [t] we show the 8-day bin cal- 
culated for the whole Fermi-LAT dataset (~ 4 years). 
As in Fig. [3] the dashed lines in Figure [7] indicate 
the corresponding to the 3cr, Aa and 5cr confidence 
level for 4 degrees of freedom. In addition to the 
three major peaks of the curve (corresponding to 
the major flares of Feb. 2009, Sept. 2010 and Apr. 
2011) we identify several episodes of enhanced gamma- 
ray emission (above a 5a pre-trial significance) that we 
call W3, W4, We and Wr in Fig. [8] These episodes have 
apparent durations in the range from 8 to 50 days (see the 
complex marked as W3 — W4), and an average fiux in the 
range = (350-500) x 10"^ phcm'^ s'^. The event 14^7 
is coincident with the very recent enhancement episode 
detected!'^ by Fermi-LAT (Ojha et al., 2012). Fig. [s] 
shows the detailed lightcurves of the Fermi-LAT wave 
episodes with the largest post-trial significance^^. 

Fig. [8] shows a remarkable episode of "wave" enhanced 
emission near MJD = 55000. The Crab was for about 
50 days above 5a from its standard gamma-ray fiux, 
with an average flux in this period F = (340 ± 6) x 
lO"'^ phcm~^ s~^. This event appears quite complex^^. 
A minor flare (that we mark as /*, at MJD ~ 54982) 
anticipates a complex and long emission that we approx- 
imate as two "waves" , W3 and W4 . The "waves" W3 and 
W4 last for ^ 15 — 20 days each, and show a post-trial 
significance (on a 8-day timescale) near or above 5a (see 
also Table 3 in the Appendix). The total episode that 
includes /* , W3 and W4 has a time duration of ^ 50 days, 
a pre-trial probability p = 2 x 10~^^ and a post-trial sig- 
nificance apost > 10. For each "wave" that we found at 
5a above the Crab average emission, we estimated the 
rise Ti and the decay time r2, the average fiux, the peak 
flux, the probability of obtaining the given x^ in the null 
hypothesis, and the post-trial signiflcance. The rise and 
decay timescales are estimated with an exponential fit. 
The results are summarized in Tab. [3l 

5. CONSTRAINTS FROM A SYNCHROTRON COOLING 
MODEL 

We can deduce important physical parameters of the 

See, e.g., |http://www. asdc.asi.it/fermi2fgl] 

^® For the Fermi-LAT data, we started witli the 2-day binned 
flux data, and then calculated the values summing over 8 days, 

= ^i—i , where Fo{i) is the i-th observed flux, Fs 

is the steady state flux, and cri is the i-th flux error. 

Due to solar panel constraints, this event was unobservable 
by AGILE. 

^® The post-trial signiflcance was calculated as in the previous 
Section. Considering four years of Fermi data, and a 8-day time 
bin, the number of trials turns out to be Nt ~ 180. 

For simplicity, we use an exponential approximation for the 
"wave" emission. Admittedly, this approximation for the episode 
of Fig. Is] centered on MJD = 55000 is not adequate. It is shown in 
Fig. [sl for illustrative purposes only. 
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Fig. 6. — Fermi-LAT 2-day bin lightcurve above 100 MeV of the Crab (pulsar plus Nebula) spanning the time period Sept. 2008/ Jul. 
2012. Results obtained with the publicly available unbinned likelihood analysis software. 
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Fig. 7. — Plot of the values (each calculated for 8-day time intervals, based on 2-day binned data) as a function of time for the 
gamma-ray Fermi-LAT data on the Crab covering the period Sept. 2008/Jul. 2012. 



enhanced 7-ray emission by adopting a synchrotron cool- 
ing model (see also Tavani et al. 2011a, Vittorini et al. 
2011, Striani et al. 2011b). Tables 1, 2, and 3 sum- 
marize the relevant information for the major "flares" 
and "waves" . We constrain the physical quantities of our 
model as described in sect. 3. We find that for flares last- 
ing 1-2 days the typical length is Z ~ (1 — 2) x 10^^ cm, the 
density constant in the range K/a = (2 — 8) x 10~^ cm~^, 
the typical Lorentz factor 7* = (2.5— 4.5)-10^, and the lo- 
cal magnetic field affecting the cooling phase in the range 
B — (1 — 2) mG. As discussed above, the total number of 
radiating particles, in case of unbeamed {5 = 1) isotropic 



emission, is iV ~ (1 - 3) • 10^^ . 

For the "wave" episodes, the typical length is I > 10^^ 
cm, the density constant is in the range K/a — (2 — 8) x 
10~^^ cm~^, the Lorentz factor 7* = (3 — 5) -10^, and the 
local magnetic field in the range B = (0.5 — 1) mG. The 
total number of particles involved is of the same order as 
in the case of flares. 

We show in Fig. |9] a schematic representation of the 
physical parameter space {B vs. 7*) for the steady emis- 
sion, waves and flares that are characterized by different 
and distinct regions. The local magnetic field for "waves" 
and "flares" is definitely amplified with respect to the 
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Fig. 8. — The most prominent wave episodes in the Fermi-LAT data with more than 5cr enhancements above the Crab average emission. 
2-day binned gamma-ray flux values (in unit of 10~® phcm"'^ above 100 MeV as a function of time. We notice that the event marlted 
as /* is intermediate between flares and waves. 



TABLE 2 

Table of the waves above Scr post-trial from the Crab average emission found in the AGILE and Fermi data. 



Name 


MJD 


Duration 


n 


T2 


Average Flux 


Peak Flux 


Pre-trial 


Post-trial 






(days) 


(days) 


(days) 


(10"**phcm-^s-i) 


(lO-^phcm-^s-^) 


p- value 


significance 


Wi 


54368-54373 


5 


2± 1 


2± 1 


440 ± 40 


670 ± 200 


4.5 X 10"** 


5.0 


W2 


54376.5-54382.5 


6 


2.5 ± 1 


2± 1 


480 ± 40 


760 ± 140 


3.0 X 10"'' 


5.5 


Wa 


54990-55008 


18 


5 ±2.5 


10 ±5 


352 ±9 


380 ± 30 


1.0 X 10"** 


4.6 


We 


55988-56000 


12 


5 ±2.5 


3.5 ± 1.5 


367 ± 12 


435 ± 35 


1.8 X 10"^^ 


6.2 


Wr 


56108-56114 


6 


3± 1.5 


3± 1.5 


431 ± 22 


450 ± 30 


1.9 X 10"'' 


5.9 



Photon fluxes are obtained for E-y > 100 MeV. 
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standard average Crab magnetic field (~ 200/iG) by a 
factor 5-10, reflecting the magnetic nature of the instabil- 
ity producing the enhanced 7-ray emission. Remarkably, 
we find that, within the monochromatic approximation 
of the enhanced particle energy distribution in agreement 
with all spectral "wave" /flare data, the Lorentz factor is 
in the range 2.5- 10^ < 7* < 5 x 10^ for both "waves" and 
"flares". The indication that both "flares" and "waves" 
have the same characteristic Lorentz factor 7* suggests a 
limitation of the acceleration process, most likely induced 
by radiation reaction. However, the monochromatic vs. 
a more extended power-law nature of the "wave" emis- 
sion requires to be tested by additional multifrequency 
data. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we addressed the issue of the duration 
and intensity of detectable enhanced gamma-ray emis- 
sion from the Crab Nebula. By considering AGILE and 
Fermi-LAT gamma-ray data above 100 MeV we find that 
the Crab produces a broad variety of enhanced emis- 
sion. We characterize this enhanced emission as short 
timescale (1-2 day) "flares" and long timescale (1 week 
or more) "waves". Given the current detection level of 
1-2 day enhancements (which can be extended to longer 
timescales of order of 1-2 weeks) , we cannot exclude that 
the Crab is producing an even broader variety (in flux 
and timescales) of enhanced gamma-ray emission. With 
the current sensitivities of 7-ray telescopes we can ex- 
plore an important but necessarily limited range of flux 
and spectral variations. It is interesting to note that 
what we called "flares" and "waves" (a somewhat arbi- 
trary division) share the same spectral properties. Given 
the current gamma-ray sensitivities, we could have de- 
tected different spectral behaviors in the energy range 
50 MeV - 10 GeV. Most likely, flares and waves are the 
product of the same class of plasma instabilities that we 
show acting on different timescales and radiation inten- 
sities. The overall detectable transient emission appears 
to be without any discernible pattern (see, e.g.. Fig. 7). 
Whether or not the instability driver of this process is 
truly stochastic in flux and timescales will be determined 
by a longer monitoring of the Crab Nebula. The tran- 
sient behavior is the topic of current intense theoretical 
investigation^'' (e.g., Bednarek & Idee 2011, Komissarov 
& Lyutikov 2011, Uzdensky, Cerutti & Begelman 2011, 
Cerutti, Uzdensky & Begelman 2012, Bykov et al. 2012, 
Sturrock & Aschwanden 2012, Clausen-Brown & Lyu- 
tikov 2012, Kohri et al. 2012, Lyubarsky 2012, Komis- 
sarov 2013, Mignone et al., 2013, Tavani 2013). The 

^'^ See the reports presented at the meeting "The Flaring Crab: 
Surprise and Impact", www.iasf — roma.inaf.it/Flaring_Crab. 



pulsar wind outflow and nebular interaction conditions 
need to be strongly modified by instabilities in the rela- 
tivistic flow and/or in the radiative properties. Plasma 
instabilities possibly related to magnetic held reconnec- 
tion in specific sites in the Nebula can be envisioned. 
However, evidence for magnetic field reconnection events 
in the Crab Nebula is elusive, and no optical or X-ray 
emission in coincidence with the gamma-ray flaring has 
been unambiguously detected to date (e.g., Weisskopf et 
al. 2012). 

Both the flaring and "wave" events can be attributed 
to a population of accelerated electrons consistent with 
a mono-chromatic or relativistic Maxwellian distribution 
of typical energy 7* ^ (2.5 — 5) • 10^. The magnetic fleld 
inducing synchrotron radiation during the decay phase 
of the waves/flares is substantially larger than in the 
steady state, as shown in Fig. [9j The range of parti- 
cle energies reproducing the wave/flare spectra is quite 
restricted with respect to the steady state, and concen- 
trated towards the maximum value of the overall distri- 
bution function. This is an important property of the 
flaring Crab acceleration mechanism whose maximum 
energy is most likely limited by radiation reaction. 

We also notice that the emitted total energies that 
can be deduced for the wave {SEj^^) and flare (SE^j) 
episodes in general satisfy the relation SEj^^ — SEjj. 
The total gamma-ray emitted energy for the wave 
episode Wi can be estimated as SE^^^i ~ lO''^ erg. Dur- 
ing the 5 days of the Wi episode, the total Crab spin- 
down energy is E^sd - 2 • lO'^^ erg. Therefore, the ob- 
served Wi efficiency for gamma-ray emission above 100 
MeV is e ~ 5 ■ 10~^. We expect the efficiency of syn- 
chrotron emission to be of order of 10% of the particle 
kinetic energy. Therefore, the energy associated with the 
wave event Wi , taken here as an example of Crab "wave" 
emission, can reach a few percent of the total spindown 
energy 

We conclude that the Crab "wave" events are highly 
significant and quite important from the energetic point 
of view. "Waves" typically imply regions larger than in 
the case of flares, and smaller average magnetic flelds. 
Their total emitted gamma-ray energy can be compara- 
ble with that associated with shorter flares. More obser- 
vations of this fascinating phenomenon are necessary to 
improve our knowledge of the flaring Crab. 
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APPENDIX 

SEARCH FOR 7-RAY "WAVE" EMISSION FROM THE CRAB NEBULA 

Table 3 shows results of our complete analysis of "wave" enhanced gamma-ray emission from the Crab Nebula. 
Seven episodes are identified in the AGILE/Fermi-LAT database with a pre-trial significance larger than 5cr. Table 
3 provides the corresponding p- values. When selected for post-trial significance (see main text), only a sub-class of 
events survives as indicated in Table 2. The detectability of the "wave" phenomenon produced by the Crab is clearly 
limited by the 7-ray sensitivity and exposure characteristics of AGILE and Fermi-LAT. Variable 7-ray emission from 
the Crab with timing and spectral characteristics different from those addressed in this paper cannot be excluded. 
Fig. 10 shows the distribution of the Fermi-LAT gamma-ray flux data of the Crab Nebula. 



TABLE 3 

Table of the waves above Sct (pre-trial) from the Crab average emission found in the AGILE and Fermi data. 



Name 


MJD 


Duration 


Tl 


T2 


Average Flux 


Peak Flux 


Pre-trial 


Post-trial 






(days) 


(days) 


(days) 


(lO^^phcm-^s-^) 


(10"Vh cm-^s-^) 


p- value 


significance 


Wi 


54368-54373 


5 


2± 1 


2± 1 


440 ± 40 


670 ± 200 


4.5 X 10"* 


5.0 




54376.5-54382.5 


6 


2± 1 


2± 1 


480 ± 40 


760 ± 140 


3.0 X 10"^ 


5.5 


r 


54980.0-54986 


6 


1 ±0.5 


2± 1 


470 ± 35 


380 ± 40 


8.0 X lO"'' 


4.2 


W3 


54990-55008 


18 


5±2 


10 ±5 


352 ±9 


380 ± 30 


1.0 X 10"* 


4.6 


W4 


55010-55025 


15 


3± 1 


6±3 


326 ± 10 


360 ± 30 


4.6 X lO"'' 


3.8 


W5 


55358-55362 


4 


2± 1 


2 ± 1 


426 ± 27 


430 ± 30 


5.6 X lO"'' 


3.7 


We 


55988-56000 


12 


5±2 


3± 1 


367 ± 12 


435 ± 35 


1.8 X 10"^^ 


6.2 


Wr 


56108-56114 


6 


3± 1 


3± 1 


431 ± 22 


450 ± 30 


1.9 X 10"^ 


5.9 



Photon fluxes are obtained for > 100 MeV. We notice that the event marked as /* is intermediate between flares 
and waves. 
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Fig. 10. — distribution of the Fermi-LAT gamma-ray fiux data of the Crab Nebula. Each value of is integrated over 8 days. 



